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The similarity between matter waves in periodic potential and solid-state physics processes has 
triggered the interest in quantum simulation using Bose-Fermi ultracold gases in optical lattices. The 
present work evidences the similarity between electrons moving under the application of oscillating 
electromagnetic fields and matter waves experiencing an optical lattice modulated by a frequency 
difference, equivalent to a spatially shaken periodic potential. We demonstrate that the tunneling 
properties of a Bose-Einstein condensate in shaken periodic potentials can be precisely controlled. 
We take additional crucial steps towards future applications of this method by proving that the 
strong shaking of the optical lattice preserves the coherence of the matter wavefunction and that 
the shaking parameters can be changed adiabatically, even in the presence of interactions. We induce 
reversibly the quantum phase transition to the Mott insulator in a driven periodic potential. 

PACS numbers: 



I. INTRODUCTION 

Since the early days of quantum mechanics toy models 
have been used to explain and teach the counterintuitive 
phenomena of the quantum world. The possibility to 
solve analytically a number of simple problems showed 
how complex processes were related to basics ideas. The 
particle in the box or the potential barrier are only two 
examples. Till the end of the last century, however, it 
was almost impossible to realize these simple toy mod- 
els because of the temperature and energy scale required 
to perform such experimental investigations. The single 
quantum object under investigation was not sufficiently 
decoupled from the external world and decoherence pro- 
cesses strongly modified the quantum evolution. The ex- 
perimental realization of quantum degenerate states with 
ultracold atomic gases allowed access to few body sys- 
tems isolated from external perturbations and operating 
at temperatures close to absolute zero. These techniques 
allowed the realization of a variety of experiments test- 
ing peculiar properties of quantum mechanics. For in- 
stance, astounding examples were achieved in the study 
of quantum tunneling in the periodic potential associated 
to optical lattices [TJ [5] , making good use of the excellent 
control on all parameters of the optical lattice defect-free 
potential. 

An important topic not extensively investigated so far 
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is that of strongly driven quantum systems, where a time 
dependent (in particular periodic) perturbation is intro- 
duced in the system in order to change its properties. 
The external driving of a parameter of the unperturbed 
system produces a modification described through an ef- 
fective potential similar to the rescaling of the mass for 
the electron motion inside a crystal. That potential con- 
tains a new variable allowing an easy understanding of 
the system evolution. In addition, the parameter tun- 
ing provides a new handle on the control of the quantum 
evolution. 

This paper presents experimental results for the mod- 
ification of quantum tunneling in a Bose-Einstein con- 
densate (BEC) using a modulated frequency/phase dif- 
ference of the optical lattice potential. By periodically 
backwards and forwards moving the spatial position of 
the optical lattice minima/maxima, the tunneling prop- 
erties of the atoms inside the lattice can be adiabatically 
changed and the atomic response is described through a 
rescaling of the tunneling rate. In addition, the ultra- 
cold atoms may be carried into novel regimes impossible 
to reach through adiabatic modifications of the Hamil- 
tonian parameters. The tunneling rescaling applies to 
a macroscopic atomic ensemble, referred to as a dressed 
matter wave in [3]. We will show how matter waves can 
be adiabatically dressed without losing the quantum co- 
herence of the ensemble, also while performing a quantum 
phase transition to the Mott insulator. Our experimental 
realization investigates simple quantum-mechanical pro- 
cesses in the driven regime and represents a quantum 
simulation of a system that is very complex to be solved 
analytically or numerically. 

The tunneling rescaling is familiar from periodically 
driven single-particle quantum systems. It occurs, among 
others, when a particle moves on a periodically forced ID 
lattice with nearest neighbor coupling such as excitons in 
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an insulating crystal [3] or electrons in a semiconductor 
superlattice [6j [7]. It also underlies the coherent de- 
struction of tunneling of a particle in a periodically forced 
double well [HI IH] , recently observed for argon atoms in 
an atomic beam [TU]. Ref. [IT] pointed out that coher- 
ent destruction of tunneling and dynamical localization 
can be interpreted as a result of destructive interference 
in repeated Landau-Zener crossings, and in this aspect 
these phenomena are similar. However, in the coher- 
ent destruction of tunneling, the initial distribution is 
frozen, while in the dynamical localization the distribu- 
tion oscillates periodically around the initial value, even 
if with a small amplitude. The tunneling rescaling is 
completely analogous to the zero order Bessel function 
rescaling of atomic g factors in the presence of oscillat- 
ing magnetic fields 1 21 examined for atoms in [13] and 
for BEC in [14] . The suppression of the Bloch band by 
the dynamical localization was observed by Raizen group 
15j for cold sodium atoms in an optical lattice when a 
weak spectroscopy probe drove transitions between the 
energy bands modified by the localization process. For a 
BEC in a periodically shaken harmonic trap, the dynamic 
splitting of the condensate, and the dynamic stabiliza- 
tion against escape from the trap were numerically stud- 
ied in [16]. A dynamical localization-like phenomenon 
occurs for the light propagation in two coupled optical 
waveguides, as predicted in [T7] and observed in [T51 119j . 
Nonlinear dynamical localization of matter wave solitons 
by means of a spatial modulation of the nonlinearity 
changing the scattering length by means of Feshbach res- 
onances and in the presence of optical lattice modulation 
was recently predicted in [2D]. 

Section 2 introduces the concept of a shaken optical 
lattice. Section 3 discusses the experimental set-up, em- 
phasizing the different techniques applied to produce the 
lattice shaking. Section 4 treats quantum mechanically 
the tunneling suppression in periodically driven systems 
for two different time dependences of the shaking force. 
Section 5 defines basic quantities characterizing ultra- 
cold gases within optical lattices. Section 6 introduces 
the Bose- Hubbard Hamiltonian for ultracold atoms in the 
presence of tunneling rescaling. Section 7 reports the ex- 
perimental results of [2D 122] on the tunneling rescaling 
in strongly driven one-dimensional quantum systems and 
on the rescaling applied to produce a Mott insulator. 



II. SHAKEN OPTICAL LATTICE 

In a ID optical lattice a standing wave is created by 
the interference of two linearly polarized traveling waves 
counter-propagating along the x axis with frequency wl 
and wavevector /cl [21 123] • The amplitude of the gen- 
erated electric field is £(r,t) — 2£Qsin(u>i 1 t)sin(ki J x). 
When the laser detuning from the atomic transition is 
large enough to neglect the excited state spontaneous 
emission decay, the atom experiences a periodically vary- 



ing conservative potential 



V 

Voi(x) = —cos(2k L x). 



(1) 



The amplitude Vq depends on the laser detuning from the 
atomic transition and on the standing wave laser inten- 
sity [53]. The periodic potential has a d^ — n/kh spac- 
ing. This potential derives from the quantum mechani- 
cal interaction between atom and optical lattice photons. 
Therefore the lattice quantities are linked to the recoil 
momentum p rcc = hk^ acquired by an atom after the 
absortion or the emission of one photon. Vq will be ex- 
pressed in units of E rec the recoil energy acquired by an 
atom having mass M following one photon exchange 



2M ' 



(2) 



We now introduce a periodic driving (often referred to 
as shaking in the following) to the system of atoms in- 
side the optical lattice. In the lattice reference frame a 
backwards and forwards motion of the periodic poten- 
tial at frequency uj along one direction is equivalent to a 
periodic force Fcos(ujt) applied to the atoms and to the 
following potential: 



V s h{t) = Kcos(uit) m\m >< m\. 



(3) 



Here K — Fd^, the modulation amplitude, is the energy 
difference between neighboring sites of the linear chain 
and |m > denotes the the localized Wannier wavefunction 
of an atom in the site in of the optical lattice [251 126] - 
The theoretical analysis of Sec. IV evidences the key role 
of the dimensionless parameter Kq defined by 



K 



K 



(4) 



III. SET-UP 

A. BEC and optical lattices 

We created BEC's of 87 Rb atoms using a hybrid ap- 
proach in which evaporative cooling was initially effected 
in a magnetic time-orbiting potential trap and subse- 
quently in a crossed dipole trap. The dipole trap was 
realized by using two intersecting Gaussian laser beams 
at 1030 nm wavelength and a power of around 1 W per 
beam focused to waists of 50 /im. After obtaining pure 
condensates of about 5x 10 4 atoms, the powers of the trap 
beams were adjusted in order to obtain condensates with 
the desired trap frequencies in the longitudinal and radial 
directions. Subsequently, the BECs held in the dipole 
trap were loaded into ID or 3D optical lattice created by 
counterpropagating Gaussian laser beams at 842 nm with 
120 /im waists and a resulting optical lattice spacing g?l= 
421 nm. The optical lattice laser beams were ramped up 
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in power in about 50 ms, this ramping time being chosen 
so as to avoid excitations of the BEC. By introducing a 
frequency difference Av between two counterpropagating 
lattice beams (using the acousto-optic modulators which 
also control the power of the beams), the optical lat- 
tice could be moved at a velocity t>l = d^Av along the 
propagation direction of those laser beams. In addition 
the optical lattice could be accelerated (or decelarated) 
with an acceleration 
F = M <2l in the rest frame of the lattice 



d^jj.Av, leading to a force 



B. Shaking set-up 

The shaking of the optical lattices was realized through 
two different schemes, described in the following. 

1) Modulated frequency difference. In this configura- 
tion the optical lattice is created by two independent laser 
beams counter-propagating along one direction. The 
beams experience separate frequency modulations in the 
passage trough two separate acousto-optic modulators 
(AOM). Each AOM is driven by a function generator 
whose frequency can be modulated both with internal 
preset functions or by using a triggered external chan- 
nel. The optical lattice shaking amplitude is related to 
the modulation amplitude Av for the frequency offset be- 
tween the two laser beams. For the case of a sinusoidal 
modulation Ai/ max sin(u>t) we have 

d 



F(t) = F max cos(u>t) = Md L — [Af max sm(o;t)] 



MdL(jjAv max cos{ut), 



leading to 



Kn = 



Hid 
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(5) 
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2) Modulated phase difference. In this case one of the 
optical lattice beams is created by retro-reflecting a laser 
beam on a mirror mounted on a piezoelectric actuator. 
The piezo-electric actuator is driven by a function genera- 
tor, and the expansion of the device is proportional to the 
voltage applied to it. For a displacement Ax max cos(ujt) 
of the mirror, we obtain 



Mu 2 Ax max rf L 



A.i 



2 w re 



(7) 



IV. DYNAMIC LOCALIZATION 

The atomic evolution in the shaken optical lattice may 
be studied by considering the localized Wannier wave- 
function | to > and the perturbations originating from 
the atomic occupation in neighboring sites [25, 26J. This 
approximation is valid when the overlap of atomic wave- 
functions introduces corrections to the localized atom 



picture, but they are not large enough to render the sin- 
gle site description irrelevant. We write for the atomic 
Hamiltonian H 

H = E ^2 \ m >< TO 

m 

- (|m >< 777 + 1| + |m+ 1 >< 777,|) . (8) 

m 

For ultracold atoms in an optical lattice of depth Vq the 
nearest-neighbor tunneling cnery J is given by |27j 



J 



-Ey, 



v rec i ^ 
l-K \ E r , 



Vn 



3/4 



exp 



E 



(9) 



The generic atomic wave function can be written as a 
superposition of the \m > localized wavefunctions 



|*(x) >=Y^C m \m > 



(10) 



The temporal evolution for the C m coefficients under the 
H Hamiltonian is given by 



ih 



dt 



EoC m + J (C m +i + C m _i) , 



(11) 



and in the following the ground state energy E will be 
supposed equal to zero. 

Dynamic localization was introduced by Dunlap and 
Kcnkrc [5 for the motion of a charged particle in the 
presence of an oscillating force, created by an electric 
field in the original formulation. It is based on exact 
calculations for the evolution on a discrete lattice and 
is valid for a generic motion within a spatial periodic 
potential. The Hamiltonian for the atomic motion on the 
linear chain with infinite sites experiencing the periodic 
potential of Eq. ^ is 

-^dyn.loc. — -J (\m >< m+ 1| + |m+ 1 ><m\) 

m 

+ Kcos (cut) \m >< m\. (12) 

in 

For the interaction with the oscillating force Fcos(ojt) 
the position operator x was assumed to be diagonal in 
the Wannier basis 128 , i.e., we have assumed 



dx < mix >< xlFxlx >< x\n >= Fd^mSm n . 



N3i 

The evolution of the wavefunction of Eq. ( 10 1 under 
-Hdyn.ioc. produces for C m the following equations: 



i — — 
dt 



J (C m+ i + C m _i) + Kcos (cot) C„ 



(14) 



The solution of these coupled equations with t = initial 
condition of atomic occupation of the m = site, i.e., 
C m (i = 0) = 6 m —o leads to [5] 

\C m (t)\ 2 = Jl (2Jt [J q 2 (Ko + Jo(K )h(t) + f 2 {t))} 1/2 ) , 

(15) 
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with the Kq parameter defined in Eq. ([6| and having in- 
troduced the first order Bessel function J m of the m— th 
order. The corresponding atomic mean-square displace- 
ment is 



y/< m 2 > 

Here 
fx(t) 
fait) 
A u (t) 

A v (t) 



V2Jt [j 2 (K Q + MK Q )f x {t) + f 2 (t))] 1/2 . 

(16) 



= 2- 



(A u + A w f 



t 2 

at (17) 
cos[Kosin(u!T)]d,T — tJo(Ko), 

Ujt 

sin[KQsin(uJT)]dT. 



The functions fi(t) and f%{t) contain the bounded func- 
tions A u and A v and therefore decay at large times. The 
expression for the |C m (t)| 2 occupation and the mean- 
square displacement are thus dominated by first term, 
unless iTm(-Ko) equals zero. At times t >> 1/oj the re- 
sulting simplified expressions of Eqs. ( 15|16 1 are 



V< m 2 > = —J e st, 
n 

where we introduced the effective tunneling rate 
Jett = JJo(Kq) 



(18a) 
(18b) 

(19) 



The applied sinusoidal force thus reduces the effective 
velocity of derealization of the initially localized atom. 
Eq. (19 1 shows that the effective tunneling energy van- 



ishes entirely whenever Kq is a root of the Jo Bessel 
function. The remarkable result is that, then, the m = 
occupation probability oscillates at frequency u> without 
decaying and that the mean-square displacement remains 
bounded. The particle is effectively localized dynamically 
by the action of the time-dependent force. 

The phenomenon of dynamic localization is shown in 
Fig. [I] in the case of time t = occupation of the 
m = site. The temporal dependence of the m = 
site is presented in (a) for K — and K — 2.405. At 
K = the decay of |Co| 2 indicates the atomic escape 
owing to the quantum tunneling. In the special case of 
the roots of the J t Bessel function the periodic recur- 
rences of \Co(t = n2Tr/uj)\ 2 , with n an integer, indicates 
that the atom returns repeatedly to the initially occupied 
site. This dynamic localization is also evident from Fig 
[TJb) where the mean displacement for the Bessel function 
roots is seen to be bounded by di,(2Tr\J\/uj) [5]. 

The modification of the tunneling rate of Eq. ( 19 ) 



appears also in Bloch periodic wavefunctions, the alter- 
native description for the quantum mechanical evolution 
within a periodic potential [55J The Bloch states 




FIG. 1: In (a) the temporal dependence of | Co j 2 plotted ver- 
sus Jt, for Ko = on the top and Kq = 2.405 on the bottom. 
For Ko equal to a root of the Bessel function, the probability 
of returning to the initially occupied site reaches unity peri- 
odically with the force oscillation. In (b) dynamic localiza- 
tion exhibited through a plot of the mean-square displacement 
V < vn? > versus Jt for different values of Ko- It is only when 
Ko is a root of Jo that the displacement remains bounded in 
time. The ratio Hcj/J = 2 in all cases. 



are energy eigenfunctions of an Hamiltonian composed 
by the kinetic energy and the periodic potential of Eq. 
([IJ and are characterized by the band index n and the 
quasi wavenumber q. For a lattice with inversion symme- 
try, with coefficients given by matrix elements between 
Wannier states located one site apart from each other, as 
in Eq. ([8]), the E n= i(q) energy dispersion takes the form 



E n=1 (q) =E + 2Jcos 



(20) 



V. ULTRACOLD ATOMS 

For a system of ultracold atoms in a ID optical lattice 
the Bose-Hubbard model, based on the tight-binding-like 
approximation [29 , considers interaction energies in a 
single site that are smaller than the gap between the 
ground state and the first excited level. While the model 
is usually based on the many body Hamiltonian with bo- 
son creation and annihilation operators at the single lat- 
tice site, for our analysis we write an Hamiltonian based 
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on the formalism of Eq. ( 12 ) 



-ffBH = -J^ (|to >< m + 1| + |m + 1 ><m\) 

m 

+ C/5> m (n m -l). (21) 

m 

Here the first sum is over nearest neighboring sites (whose 
number is determined by the one, two or three dimension 
geometry), and n m is the number operator that counts 
the number of atoms at the m-th site. Notice that this 
approach cannot be applied to describe the Mott insula- 
tor dynamic for which the second quantization approach 
is required. The dynamics of cold atoms in deep lattices 
is expressed in the Bose-Hubbard model by using only 
two parameters: the tunneling energy J and the on-site 
energy U of two atoms occupying the same lattice site, 
i.e. the energy required for the presence of more than 
one particle per site. By using simple approximations on 
the wavefunction [27] , the dependence of U on the lattice 
depth Vo may be calculated. In the ID experimental re- 
alization it is necessary to consider the orthogonal direc- 
tions along which the lattice is not present. The energy 
contribution to the system can be evaluated by writing 
the wavefunction for the transverse directions as in the 
weak harmonic potential that in our case is due to the 
dipolar trap confinement of the cloud. The ID parameter 
U becomes 



U : 



3. 



1/4 



'y^z 



E IC 



(22) 



where a s is the atomic scattering length and w yjZ are the 
harmonic trap frequencies in the transverse directions. In 
3D optical lattice with compression applied by the light 
in all directions, U becomes 



U = —=ka s E^ 



V 



3/4 



(23) 



problem is to use the Floquet theory 30J , or the dressed 
atom theory |12j . The solution of the Schrodinger equa- 
tion with the Hamiltonian of Eq. ( 24 1 has solutions of 
the form 



|Vn(t) = |«n(*) > 



(25) 



where the so-called Floquet mode \u n (t) >= \u n (t + T) > 
is again periodic in time with period T. The energy e n is 
called a quasienergy because of the formal analogy with 
the quasimomcntum in the Bloch problem in a spatially 
periodic Hamiltonian. By substituting the Floquet solu- 
tion of Eq. ( 25 1 into the Schrodinger equation one ar- 
rives at an eigenvalue problem [3T]. Notice that the 
\ u !n,o)(t) > is a solution of that problem with eigenvalue 
e n , while then |u( n ,m)(*) > = l u (n,o)(*) > exp(imujt) is 
also a solution with eigenvalue e n + mhw, where m is any 
positive or negative integer. At times t = to + sT with 
integer s, they all coincide, apart from a phase factor. 
Hence the Floquet spectrum repeats itself periodically 
on the energy axis. Each quasi-energy band of width huj 
contains one representative, labeled by m, of the class of 
eigenvalues belonging to the Floquet state labeled by n. 

Ref. [32j [33] demonstrated that the presence of a driv- 
ing force corresponds to the following rcnormalization of 
the Bose-Hubbard Hamiltonian: 

H xen = -J e ffV] (|m >< m + 1| + \m + 1 >< m\) 

m 

+ L/^n m (n m -l). (26) 

m 

where the tunneling rate J is substituted by the effective 



tunneling rate J c g- defined in Eq. ( 19 ) and the interaction 
energy U is not modified by the shaking. 

This renormalized description applies also to the case 
when an additional static force is applied to the atoms, 
(theory in [34 and experimental realizations in |35[ 136]). 



VI. BOSONS IN PERIODICALLY DRIVEN 
OPTICAL POTENTIAL 



VII. EXPERIMENTAL RESULTS 



We now add a periodic driving to the Bose-Hubbard 
model, with the atomic motion described by the Hamil- 
tonian 

H = Hsa. + Kcos(u>t) to | to >< m\ 

m 

= - (\m >< m + 1| + |777 + 1 >< 771 1) 

m 

n m - 1) 

m 

+ Kcos(ujt)^^m\m><m\. (24) 
m 

The full Hamiltonian is now periodic in time with pe- 
riod T = 2ir/u} and a good strategy for this theoretical 



In a preliminary experiment without shaking (Kq = 0), 
we verified that, for our expansion times, the growth in 
the condensate width a along the ID lattice direction 
was to a good approximation linear and for a/di, in very 
good agreement with the theoretical dependence of Eq. 
(18b I with the tunneling given by Eq. (J9j> . Experimen- 



tal results for a versus the expansion time t at a given 
depth of the optical lattice are plotted in the bottom part 
of Fig. 2(a), and an image of the expanded condensate 
cloud at t = 150 ms is reported in Fig. 2(b). This enabled 
us to confirm that da / dt measured at a fixed time was di- 
rectly related to J and, in a shaken lattice, to | J c s (Kq) |. 
The expansion of the condensate width versus time at 
modulation parameter Kq = 2.4 corresponding to the 
first zero of the Bessel function, is plotted in Fig. 2(a). 
The corresponding condensate image is in Fig. 2(c). The 
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dynamical localization leads to a blocking of the conden- 
sate expansion. By varying the optical lattice depth Vq 
and the shaking frequency u> we verified that the univer- 
sal behavior of \J c s/J\ was in very good agreement with 
the zero-order Bessel function rescaling of Eq. (19) for 
Kq up to 6. 




(a) 



50 100 150 

Expansion time (ms) 
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FIG. 2: Results for the free expansion of the condensate, mea- 
sured in situ. In (a) plot of the condensate width versus versus 
expansion time at Ko = (open square) and Ko = 2.4 (closed 
square) at Vo = 6-E rC c- The straight lines fitted through the 
data allowed us to derive a reduction of the tunneling by a 
factor around 25. In (b) and (c) images at expansion time 
t =ms in a Vo = E lec lattice, for different values of the shak- 
ing parameter, Ko = and Ko = 2.4 respectively. 

Within the range 2.4 < Kq < 5.5 the zero-order Bessel 
function changes its sign, and also the tunneling rate. We 
verified this sign change by monitoring the phase coher- 
ence of the BEC in the shaken lattice, which was made 
visible by switching off the dipole trap and lattice beams 
and letting the BEC fall under gravity for 20 ms. The re- 
sulting spatial interference pattern is a series of regularly 
spaced peaks at 2n x p ICC with n integer (positive or neg- 
ative), corresponding to the various diffraction orders of 
the propagating matter waves. A standard interference 
pattern, as corresponding to the J value of the tunneling 
rate at Kq = 0, is presented in Fig. 3(a) with the maxi- 
mum at p — and two secondary maxima at p — ±2p rcc . 
Such a pattern may be interpreted as the multiple source 
interference of Bloch waves with quasi-momentum q = 
extending over the whole optical lattice. In the region 
between the first two zeros of the Bessel function, where 
J off < 0, we found a typical interference pattern as in 
Fig. 3(c). That interference pattern is produced by a 
staggered Bloch wavefunction with q — ±k^, at the edge 
of the Brillouin zone [37 . This different condensate wave- 
function is produced by the inversion of the curvature of 
the (quasi) energy band at the center of the Brillouin zone 



when the effective tunneling parameter is negative, with 
energy minimum at the Brillouin zone edge, as plotted 
in Fig. 3(d). 



(b) 




(d) 
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FIG. 3: In (a) and (c) interference of dressed matter waves 
released from an optical lattice with Vo = 9E IBC and oj/2-k = 
3 kHz. The interference was measured by switching off the 
dipole trap and lattice beams and letting the BEC fall under 
gravity for 20 ms. In (a), interference at Ko ~ 1.5, corre- 
sponding to J c fr/ J = 0.51; in (c), Ko = 3 corresponding to 
Jes/J = —0.26. In (b) and (d) calculated ground energy 
band structure E\(q) versus the quasimomentum q for the 
corresponding values of Vo and Ko- 

In [32] we realized the coherent control of the dressed 
matter waves. Thus we changed adiabatically and re- 
versibly the quantum state of ultracold bosons in driven 
optical lattices between a superfluid and a Mott insulator 
(MI) by varying the amplitude Kq of the shaking. Within 
the Bose-Hubbard model described by the J and U pa- 
rameters, if U/J << 1 the tunneling dominates and the 
atoms are delocalized over the optical latice. Instead for 
U I J > > 1 the interaction term leads to a loss of phase co- 
herence through the formation of number squeezed states 
with increased quantum phase fluctuations. At a criti- 
cal value of U/J the system undergoes a quantum phase 
transition to a MI state [38]. Using optical lattices one 
can tune U/J by changing the lattice depth, which af- 
fects both U and J through the width of the on-site wave 
functions. In addition at fixed lattice depth U can be 
increased making use of a Feshbach resonance. Alterna- 
tively we have suppressed J by periodically shaking the 
lattice and verified that the MI state is reached for a 
critical value of U / J e s- 

In order to realize the driving-induced superfluid-MI 
transition, we first loaded a BEC into a 3D lattice with 
Vo = H-E-roc using an exponential ramp of 150 ms dura- 
tion and then linearly increased Kq from to Kq = 1.62 
in 4 ms, as schematized in the lower part of Fig. 4. 
While in an undriven lattice with at HE rcc lattice depth 
the BEC is superfluid with {7/6 J — 3.5, as shown by the 
interference pattern in the top lef of Fig. 4, for the driven 
lattice at Kq = 1.62, U/6J e g — 7.9, i.e., larger than the 
MI transition critical value. The distinct loss of phase 
coherence observed into the interference pattern on the 
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top right of Fig. 4 confirms that the system was in the 
Mott insulating phase, as in [38]. When Kq was ramped 
back to 0, the interference pattern reappeared, proving 
that the transition was induced adiabatically and that 
the system was not excited by the driving. 




time (ms) 



FIG. 4: Shaking-induced Mott insulator transition, (a) In a 
3D lattice of constant depth (Vb = lli?rec, as in the right 
scale) and driven at u>/2ir = 6 kHz, Ko was ramped from to 
Ko = 1.62 in 4 ms and back, as schematized by the lines with 
the left scale. The interference patterns observed at Ko — 0, 
at Ko = 1.62 and back at Ko — are included. 

In addition we compared the visibility of the interfer- 
ence pattern by inducing the MI transition in two dif- 
ferent ways: (a) by increasing Vb in an undriven lattice 
as in[3S] and (b) by varying K for constant Vq. The 
dependence of the visibility on U/6J c s is the same for 
methods (a) and (b), strongly indicating that the same 
many-body state is reached by the conventional approach 



of increasing U or by our approach of keeping constant 
and reducing J through the dynamical localization. 



VIII. CONCLUSIONS 

Our results confirm and extend the role of cold atoms 
in optical lattices as versatile quantum simulators and 
open new avenues for the quantum control of cold atoms, 
thus establishing a link to coherent control in other sys- 
tems such as Cooper pairs in Josephson qubits. The ex- 
plored scenario is not intended as a look at the common 
superfluid-insulator transition from a different angle, but 
aims at obtaining genuinely new, nontrivial information 
on condensate dynamics. The control demonstrated here 
can be straightforwardly extended to more than one driv- 
ing frequency and to more complicated lattice geometries 
such as superlattices [33]. In addition ref. @U] suggested 
that a control on the phase of the driving field may be 
used to produce and maintain a coherent atomic current. 
The application of dynamical localization to a triangu- 
lar lattice Bose-Hubbard model allows the modeling of 
geometrically frustated antiferromagnetim [4d] . 
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